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We propose and analyze a new approach to the coherent control and manipulation of quantum 
degrees of freedom in disordered, interacting systems in the many-body localized phase. Our ap¬ 
proach leverages a number of unique features of many-body localization: a lack of thermalization, 
a locally gapped spectrum, and slow dephasing. Using the technique of quantum phase estimation, 
we demonstrate a protocol that enables the local preparation of a many-body system into an ef¬ 
fective eigenstate. This leads to the ability to encode information and control interactions without 
full microscopic knowledge of the underlying Hamiltonian. Finally, we analyze the effects of weak 
coupling to an external bath and provide an estimate for the fidelity of our protocol. 

PACS numbers: 71.55.Jv, 05.30.Rt, 64.70.P-, 72.20.Ee, 71.23.An 


The coherent control and manipulation of a complex 
quantum system is one of the central challenges of mod¬ 
ern physics. The majority of ongoing research focuses 
on building up this complexity starting from individual, 
isolated qubits [1-7]. In contrast, the opposite approach, 
where one seeks the coherent manipulation of a strongly 
interacting system, especially subject to disorder, is gen¬ 
erally thought to be intractable. The main difficulty 
is the exponentially growing Hilbert space with a typ¬ 
ical many-body eigenstate strongly coupled to a dense 
set of other states. Even when the structure of relevant 
eigenstates is known, their many-body character gener¬ 
ally makes them difficult to manipulate using most ex¬ 
perimental controls. 

In this Letter, we explore an alternate approach to¬ 
ward the coherent control of many-body systems. Our 
approach follows a new paradigm introduced by recent 
studies of many-body localization (MBL) [8-33]. In the 
presence of strong disorder, the many-body eigenstates of 
an isolated, interacting system can be localized in Fock 
space [13]. These MBL eigenstates exhibit a discrete lo¬ 
cal spectral response, suggesting that one can coherently 
manipulate and store quantum information (using exter¬ 
nal controls with both finite spatial and spectral resolu¬ 
tion) [21, 22, 24]. This is remarkable because it implies 
that, in the presence of sufficiently strong disorder, local 
quantum bits naturally emerge from a strongly interact¬ 
ing many-body system [22, 25]. 

Further, in the MBL phase, these emergent qubits in¬ 
teract with one another via diagonal (i.e., dephasing) 
interactions; for two qubits separated by a distance d, 
this interaction time scales as ^ rexp(d/^), where r is 
a characteristic timescale and ^ is the localization length 
[16, 20]. In this respect, the MBL phase differs crucially 
from, e.g., a conventional Anderson insulator, in which 
such interactions are absent [8, 14, 34]. This weak inter¬ 
action brings about two implications for quantum con¬ 
trol: (i) the coherence time of a single MBL qubit can 
be extremely long and (ii) dephasing can be leveraged 



FIG. I. (a) Schematic energy spectrum of the disordered 
spin-1 model, (b) A discrete local spectrum with excitations 
exhibiting different localization lengths, (c) External elec¬ 
tromagnetic driving (yellow beam) defines the location of an 
MBL qubit. 

as a way to coherently couple spatially separated MBL 
qubits [16]. Unfortunately, in general, the parametric 
dependence of both the coherence and dephasing times 
are identical. To make use of many-body localization 
as a platform for quantum control, one must devise a 
method to implement interactions (e.g. two-qubit quan¬ 
tum gates) on much shorter timescales than the decoher¬ 
ence. 

Our approach relies upon three key ingredients. First, 
the local spectral gap of the MBL phase enables one to 
encode and manipulate quantum information using con¬ 
trols of both finite spatial and spectral resolution. Sec¬ 
ond, we demonstrate the ability to enhance interactions 
between MBL qubits by exploiting variations in the local¬ 
ization length of different eigenstates. Finally, we utilize 
quantum phase estimation and projective measurements 
to locally prepare the system in a many-body eigenstate 
ensuring that the interactions between spatially sepa¬ 
rated MBL qubits is coherent. In contrast to previous 
work [28], our approach uses the intrinsic dephasing in¬ 
teractions of the MBL phase as a resource to purify the 
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system. The combination of these ingredients can in prin¬ 
ciple, form the basis for a many-body localized quantum 
information processor. We find that such an approach 
can achieve a fidelity that asymptotically approaches the 
conventional method using isolated qubits. Perhaps more 
significantly, the local encoding of information provides 
new routes for studying and quantifying dynamics in the 
MBL phase. For example, one can imagine directly prob¬ 
ing the decay of both the classical and quantum informa¬ 
tion, extracting the characteristic decay length for each 
[35], thereby gaining new insights into the nature of MBL 
systems. 

Spin-1 model —To illustrate the key idea of our ap¬ 
proach, we start with a simple model Hamiltonian that 
clearly exhibits two separate localization lengths. We 
will utilize excitations with a short localization length 
(“qubit” excitation) to encode information; meanwhile, 
to implement gates, we envision off-resonantly coupling 
such qubit excitations to a more delocalized excitation 
(“bus” excitation), thereby enhancing interactions on de¬ 
mand [4, 36]. 

Consider a chain of spin-1 particles with Hamiltonian, 

N 

H = Y, hiSt + Ai {Sff + 2tSi ■ Si+i (1) 

i=l 

where N is the total number of sites (Fig. la). Disorder 
exists in both the magnetic field strength hi and the sin¬ 
gle spin anisotropy A^, which are randomly drawn from 
uniform distributions [ho —h, ho+h] and [Aq —A, Aq+A], 
respectively. As the Hamiltonian conserves total z- 
magnetization, define a logical ground 

state \G) as the eigenstate with all spins in \Sf = — 1). 
We focus on the regime in which the field disorder is much 
weaker than the anisotropy disorder, with a hierarchy, 

t < hAo/A < A < Ao < ho. (2) 

In this limit, there are two types of well-defined single¬ 
particle excitations. The first (qubit-type) consists of a 
single spin flipped into the \Sf = 0) state. This excitation 
moves in a disorder potential of strength A, analogous 
to Anderson localization [8]. As such, all single-particle 
states are localized and in the limit t ^ A, the typi¬ 
cal localization length of an excitation is given by ^ 
log A/t. The second single-particle excitation (bus-type) 
consists of a single spin flipped into the \Sf = +1) state; 
it can be regarded as a bound state of two qubit-type 
excitations, with binding energy ^ Aq (Fig. la). A bus 
excitation propagates via a second-order hopping pro¬ 
cess of amplitude ^ t^/(2Ao). It, too, can be regarded 
as a localized single-particle excitation, with localization 
length, ^ 2 “^ log (hAo/e) = - log [tA/(hAo)]. Un- 

der the hierarchy of equation (2), ^2 ^ (Fig- lb). 

We now turn to the case of multiple excitations. For 
sufficiently small t or sufficiently low energy densities. 


the many-body eigenstates of this system are expected 
to exhibit many-body localization. The effective inter¬ 
action between two excitations in the MBL phase falls 
off exponentially with distance; in the limit of low en¬ 
ergy density, it is simply given by the overlap of their 
wavefunctions. Thus, for two qubit-type excitations a 
distance d from one another, the interaction strength is 

<5/~ texp(-d/^i), (3) 

while for two bus-type excitations, the interactions scale 
as ^ texp(—(i/^ 2 )- 

Single qubit eneoding and manipulation —In order to 
encode Nq qubits, we partition a given disordered sys¬ 
tem of size L into Nq = L/d smaller segments of length 
d > ^2 ^ ^ 1 - We select a single qubit-type excitation 
centered near the middle of each segment, and use it to 
encode the state of the MBL qubit (Fig. Ic). An arbi¬ 
trary superposition of qubit states will dephase on a time 
scale 1/d/, and all gate operations must be fast compared 
to this timescale [16, 20]. 

The coherent manipulation of a single qubit state is 
feasible with local external control of finite spatial £ and 
spectral 6uj resolution [24]. After initializing the sys¬ 
tem to the logical ground state |G), one applies a time- 
dependent electromagnetic field of frequency uj (Fig. Ic); 
to within exponentially small corrections, this transverse 
field couples \G) with only a finite number ^ max(^, ^ 1 ) 
of localized eigenstates, as depicted in Fig. lb. 

Thus, one can induce a high-fidelity transition between 
\G) and a specific eigenstate by weakly driving the sys¬ 
tem with (j tuned to the particular many-body transition 
ujl- The requirements for the spatial and spectral reso¬ 
lution are: (i) ^ < d and (ii) Sou < A/max(^,^i). In 
this encoding scheme, information readout is performed 
by directly measuring the local spin polarization. In the 
limit ^ a, this corresponds to the measurement of a 
single spin. When > a, the total polarization in the 
region determines the qubit state. 

Two-qubit gates —To perform universal quantum con¬ 
trol, we only need to demonstrate a controlled phase gate. 
As previously mentioned, such a gate can in principle, 
be achieved via the dephasing interactions of Eq. (3). 
However, since decoherence occurs at the same rate, one 
must effectively enhance this interaction strength. Since 
the dephasing between bus excitations occurs at a much 
faster rate te~^l^^ ^ the effective interaction strength can 
be enhanced by several orders of magnitude upon dress¬ 
ing [37] a qubit-type excitation with a bus-type excita¬ 
tion. This enhancement ratio is given by. 



where Q is the Rabi frequency and 6 is the detuning of the 
external driving field. A numerical demonstration of this 
enhancement in the simplified spin-1 model is provided 
in the Supplemental Material [38]. 
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FIG. 2. Numerical study of localization lengths in the low spin 
density, weakly localized regime of Eq. (6). (a) Histogram 

of localization lengths in units of lattice spacing. Since the 
model conserves total , one can choose to work in manifolds 
with a fixed density of excitations. Here, we consider two spin- 
flips. (b) Decay of the interaction strength as a function of 
distance for both qubit and bus excitations. 


Generic MBL system —We now generalize our protocol 
from the spin-1 model to any MBL system with a fully 
localized spectrum [8-27]. Such systems can be repre¬ 
sented in terms of conserved local spin-1/2 degrees of 
freedom termed “1-bits” [16, 20]: 


HmU = ' ' ' ( 5 ) 

i i,j 

where the are spin-1/2 operators describing localized 
excitations of size hi is the energy of each excitation, 
Jij is the two-body interaction shift, etc. An 1-bit over¬ 
laps with ^ microscopic degrees of freedom, and can 
therefore be individually addressed (up to exponentially 
small infidelities) using a drive with spectral resolution 
greater than t/2^. In addition, the 1-bits have effective 
diagonal interactions that also fall off exponentially, with 
a localization length ^ that varies from site to site and 
(unlike the 1-bit size from eigenstate to eigenstate [20]. 
All of the qualitative arguments from the spin-1 model 
extend to the generic case, provided (a) there are excita¬ 
tions that overlap spatially but have different localization 
lengths and (b) one can initially prepare the system in 
a particular eigenstate. 

To understand the first of these issues, we numerically 
study the variation of localization lengths for the Hamil¬ 
tonian, 


H = Y,2tS^■S^+^ + h^St (6) 

i 

where hi is randomly drawn from a uniform distribution 
[—A, A]. This model is “generic” in that we have not 
engineered its localization lengths to have any particu¬ 
lar features. In Fig. 2a, the histogram of localization 
lengths ^ is shown for low excitation density (up to two 
spin-flips) with weak disorder A/t = 3 [38]. The distri¬ 
bution exhibits heavy tails, evincing the fact that there 


is significant probability to find ^ longer than its average 
value. It turns out that variations always exist (even at 
high excitation densities) and are largest near the local¬ 
ization transition (see [38] for a quantitative analysis). 

We note that our protocol does not require a para¬ 
metric separation between long and short localization 
lengths; rather, one can optimize d to achieve as large 
an interaction enhancement as desired. To quantify the 
enhancement of the interaction strength, we numerically 
extract the localization length of typical qubit and bus 
excitations. We diagonalize systems up to length L = 50 
with at most two excitations and perform 100 disorder 
averages. For each disorder realization, we apriori choose 
approximate positions for the left and right MBL qubits 
(at L and R). Starting from the polarized state |G), 
we identify computational qubits with significant local 
support around L, R. These basis states are chosen as 
energy eigenstates that have significant overlap with local 
excitation operators ^l(r) 

In particular, we select states (0^,0^,0^^), such 
that |(0^|E^|G)|M(0^|E^|G)P, |(0^^|E^E^|G)| > 

10“^; moreover, to ensure that the MBL 
qubits are self-consistently defined, we require 

> 10 - 2 . 

This enables us to directly estimate the effective inter¬ 
action between the MBL qubits from energy eigenvalues, 
with 51 = Elr — El — Er + Eq- As the qubit (bus)-type 
excitation, we select basis states, (0^, 0^, 0^^) that ex¬ 
hibit the smallest (largest) interaction strength [38]. As 
depicted in Fig. 2b, bus-type excitations interact signif¬ 
icantly more strongly than qubit-type excitations, en¬ 
abling multiple quantum gates to be performed within 
a decoherence window. While this simplified simulation 
does not account for the dynamical “dressing” of the 
qubits, the nature of the qualitative enhancement is the 
same: it comes from intrinsic variations of the many- 
body localization length (see [38] for full numerical sim¬ 
ulations). 

Initial state preparation —To enable fully coherent evo¬ 
lution, the interactions between MBL qubits must be co¬ 
herent. This is indeed the case when the system is pre¬ 
pared in a fiducial many-body eigenstate; however, such 
a preparation is very difficult for a generic (high tempera¬ 
ture) system. Since the MBL phase has a simple descrip¬ 
tion in terms of local conserved quantities, each many- 
body eigenstate can be labeled by specifying the values 
(±1) of all r spins. The question is then: how can one effi¬ 
ciently prepare the system into a desired r-spin superpo¬ 
sition, starting from an arbitrary many-body state, with 
only local control over small sub-regions A and B. Let us 
label the two excitations L and R (in each region) to be 
our MBL qubits. The effective interaction between these 
qubits takes the form 6i = Jlr + ^LUR^k + • • • 

for each r-spin configuration [20]. Thus, starting from a 
generic many-body state, the effective interaction varies 
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from eigenstate to eigenstate and cannot enable a coher¬ 
ent quantum gate. 

Interestingly, the interaction itself can be used as a 
resource to purify the entire system of r-spins; in par¬ 
ticular, by observing the interaction strength between 
the MBL qubits, one can effectively perform a quantum 
non-demolition measurement of the r-spin configuration. 
Such an observation can be done via a modified spin-echo 
protocol [34], which projects the r spins onto a set of con¬ 
figurations that have the same up to the precision of 
the measurement. 

Specifically, using an adaptive phase estimation algo¬ 
rithm [2, 28, 39], one can repetitively measure the inter¬ 
action strength, 5i = 27rt 5c,2“^ (in binary), to a 

precision set by its smallest significant digit 5 m- To mea¬ 
sure each digit, 5/^, requires a time, T/^ = ^, yielding a 
total measurement time, 

M 

Ttot = J2^k = 2TM{l-2-^). (7) 

k=l 

A few observations are in order: (1) this procedure is 
extremely efficient since the number of measurements M 
scales logarithmically with the desired precision and (2) 
the total measurement time Ttot is also the time-scale 
over which the MBL qubits can now be expected to inter¬ 
act coherently. This result implies that one can perform 
coherent quantum manipulations in the infinite temper¬ 
ature MBL phase with a preparation overhead scaling 
only linearly in time. 

Imperfections —In what follows, we analyze a variety of 
realistic imperfections and provide a quantitive estimate 
for achievable fidelities in a number of experimental sys¬ 
tems. In particular, we will consider the effect of finite 
spatial(^) and spectral((5cj) resolution, leading to: (1) im¬ 
perfect initialization, (2) off-resonant excitation, and (3) 
population loss into nearby modes. 

For single qubit gates, the fidelity is given by, 

where T 2 characterizes the extrinsic decoherence time 
and £ represents the local spectral gap. The lack of in¬ 
dividual addressability manifests as a ratio of the spatial 
extent to the qubit localization length while off- 

resonant excitations induce an error i/{T2e?. 

Turning to the fidelity for two qubit gates, we note that 
finite resolutions (spatial and spectral) bring about two 
consequences, namely, an optimal choice of qubit sep¬ 
aration dopt and a renormalized bus localization length 
^2 ^ [38]. The former arises when the decoherence 

rate begins to dominate the bare interaction strength, 
while the latter occurs for line-widths larger than the 
local spectral gap. The physics of this latter case is anal¬ 
ogous to coupling a single bound state (qubit excitation) 


with a multi-particle continuum (bus excitation) of ef¬ 
fective mass m* 1/ta^, and leads to a renormalized 
excitation size ~ a\/tT 2 [38, 40-42]. Combining these 
two effects gives an overall fidelity, 

J -2 - 1 - (9) 

where F P/tTfa^ + I/T 2 is an effective decoherence 
rate estimated from Fermi’s Golden rule and ^ 2 ^^ — 
^ 2(^0 spectrally (un)-resolved bus excitations. 

Interestingly, these bounds are consistent with tradi¬ 
tional quantum information processing schemes based 
upon isolated qubits. In particular, for ^ 1 ,^ ^ a and 

^2 ^ 00 , one discovers 1 — (tT 2 ) , This 

corresponds to two tightly localized qubits (e.g. atoms or 
bound states) interacting via a band of delocalized states 
(e.g. phonons or photons). 

Fidelity Estimates —Our protocol applies most read¬ 
ily to quantum optical systems with local addressing. 
A number of such platforms are promising candidates 
for realizing many-body localization, including ultracold 
atoms, dipolar molecules, superconducting qubits, and 
solid-state spins [43-46]. In the case of ultracold atoms, 
a direct implementation of a spin system is feasible via 
multi-component Fermi- or Bose-Hubbard models. From 
recent experiments [47-49], the spatial resolution ^ ^ a, 
the typical superexchange interaction strength t ^ lOHz, 
and the coherence limited by particle loss T 2 ^ 3s are 
feasible, yielding an overall fidelity IF 2 ~ 0.92. Recent 
progress towards the engineering of large superconduct¬ 
ing flux-qubit arrays is particularly intriguing [50, 51]; 
disorder naturally arises from the fabrication process and 
full tomography of the couplings within the system is 
daunting. Thus, the ability to define MBL qubits in a 
modular fashion is particularly applicable. With recent 
coherence times [52] > lO/rs, typical interaction strengths 
^ IGHz, and individual flux-qubit control, one finds a 
fidelity J ^2 ~ 0.99. In the case of molecules and solid- 
state spin impurities, the interactions are long-range and 
the dominant disorder arises from random bonds. For ef¬ 
fectively short-range power-laws, many-body localization 
persists and the main idea of this work is still applicable 
[26, 29, 30, 43]. 

In summary, we have introduced a scheme for the co¬ 
herent control of local degrees of freedom in the many- 
body localized phase. Our approach enables encoding 
quantum information as well as to perform quantum logic 
between separated MBL qubits. This suggests that in 
certain cases, strongly disordered, interacting systems 
may be a resource for quantum information applica¬ 
tions. The ability to efficiently prepare (high temper¬ 
ature) many-body eigenstates via local spectroscopy also 
opens the door to studying coherent dynamics in the 
MBL phase. By probing the decay of both classical and 
quantum information, it may be possible to characterize 
many-body localized states and their dynamics. 
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Numerical methods and details —Here, we provide full numerics for our protocol in both the spin-1 model in 
Eq. (1) and the generic Heisenberg spin-1/2 model with random fields (Eq. 9). The general outline of the demonstration 
is the following. Eirst, for a given realization of a disordered Hamiltonian we numerically identify two qubit- and 
bus-type excitations, located at positions L and R. Second, we compute the effective interaction strength between 
qubits either with or without off-resonant coupling to the bus excitations. Third, we repeat this process for various 
distances separating L and R and average the interaction strength over 500 realizations. Below, we provide the details 
of this process for both spin-1 and spin-1/2 models, always setting the lattice spacing a = 1. 

Eor the spin-1 model, we exactly diagonalize a system of 14 spins (Eq. 1) with Aq = 25, A = 15, ho = h = 0.02, 
and t = 1. The Hamiltonian conserves total z-magnetization, Therefore, we only need to diagonalize a few 

relevant symmetry sectors of the Hilbert space, namely, systems with up to four spin excitations starting from the 
spin polarized state \G). A single qubit excited state at position i is identified as the energy eigenstate with the 
largest overlap with Sf \G) where Sf is the spin flip operator at position i. We denote the eigenstate with a single 
excitation at position L as |0^). In this simulation we choose L = 3. Similarly, we identify eigenstates with single 
qubit excitations at various positions i? = 8, 9, 10, 11, and 12th site and denote them as |0^). The eigenstates with 
two qubit excitations and bus excitations at L{ot R) are also identified in a similar way from \ <p^) and 

|0^)(or respectively. The off-resonant driving on the left qubit has been implemented by transforming 

the Hamiltonian into the rotating frame with a rotating wave approximation. The time-dependent driving on the 
right qubit has been numerically integrated over one period, and then the long time dynamics is studied in the Eloquet 
basis. As discussed in the main text, the effective interaction strength between the MBL qubits is directly calculated 
from energy eigenvalues. 


5i = Elu — El — Eji + Eq^ (SI) 

This is also the case for the driven hamiltonian in the Eloquet basis, where energies are defined via logarithms of 
unitary evolution. The average interaction strength as a function of distance is shown in Eig. Sla. We confirm the 
exponential decay of the bare interaction strength 5i with a length scale 0.53 (blue solid line), which agrees with 

our simple theoretical estimate (log y) 0.5. When the qubits are driven, the interaction is enhanced by several 
orders of magnitude (red dash-dot line and green dotted line). Moreover, we observe that the enhanced interaction 
strength scales quadratically in ^2/(i (see Eig. Sib), confirming the off-resonant dressing picture presented in Eq. (7) 
of the main text. 

We use a similar technique to simulate a system of 16 spins in the Heisenberg model, with t = 1 and A = 7.5. 
Again, the simulation supports up to four excitations from the polarized state \G). Unlike our spin-1 model, the 
qubit excitations are identified by their energy eigenvalues; single spin excitations with the highest and the lowest 
energies are used as |(/)^) and depending on their relative positions. A single excitation with either high or 

low energy is expected to have a short localization length. Therefore, coupling to any subsequent excitations will, 
on average, enhance the effective interaction. Note that our criteria of identifying qubit- and bus-excitations can be 



FIG. SI. Effective interaction strength, (a) Interaction strength decays more weakly when either one or both qubits are 
off-resonantly coupled to the bus excitations, (b) The enhanced interaction strength is quadratic in the ratio of Rabi frequency 
to detuning Q/(5 as expected for off-resonant dressing, (c) Effective interaction strength decay for the random field Heisenberg 
model. 
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FIG. S2. Variation in the localization lengths. Left: the phase diagram of the random held Heisenberg spin-1/2 model as a 
function of spin density and disorder. Middle: histograms of the inverse localization length for parameters shown in the left 
panel. Right: mean and standard deviation of the inverse localization length in the strong disorder regime. 


further optimized to achieve a stronger interaction enhancement. We estimate the positions of qubit excitations Pl(r) 
by maximizing | \G) |, and then identify bus excitations from Sp^_^i and |0^), i-e. nearest 

subsequent excitations. The effective interaction strength as a function of distance is shown in Fig. Sic. Again, we 
conhrm that the enhancement is exponentially large at long distances. We emphasize that this model does not have 
‘engineered’ localization lengths as in the spin-1 model. 

The typical localization lengths of qubit and bus excitations can be obtained from the slopes in Fig. Sic, yielding 
^ 0.34 (blue solid line) and ^2 ^ 0.50 (red dash-dot line). These values are consistent with the probability 
distribution of localization lengths studied in the next section. 


Variation of localization lengths —Here, we provide numerics showing that a variation in localization lengths 
always exists in the localized phase of the random field Heisenberg model (this variation is of course most pronounced 
near the transition). 

First, to work in the localized phase, we estimate the criticial disorder A/t for the many-body localization phase 
transition for various spin excitation densities (relative to the polarized state |G)). We follow the method used in Ref. 
[1], and calculate the fraction of dynamical polarization at infinite temperature. The result is summarized in the left 
panel of Fig. S2. We confirm that the transition occurs at A/t ^ 5 for high densities, consistent with previous results 
for small system sizes [1]. At low densities, the critical disorder is smaller, approaching the limit of the Anderson 
localization criteria in ID A/t ^ 0. 

We study the variation of localization lengths under various conditions. Here, we define the localization length of 
excitations with respect to a particular a many-body eigenstate in the following way. First, we exactly diagonalize 
a system of 12 spins, choosing an arbitrary energy eigenstate with total spin up density + |))/^ 

define it as our logical ground state \Q) = |00). Then, we choose two “probing positions” Lp = 2 or 4 and Rp = 9 
or 11 and identify all possible computational basis states {viz. |01), |10), and |11)) as energy eigenstates that have 
large enough overlap with |00), |00), and |00); we also impose that these states satisfy the consistency 

condition | (111 |01) |, | (11| |10) | > c with a threshold c = 0.01. The relative differences in the energy eigenvalues 

of these states defines the interaction strength and the localization length ^ by 


d/e 



= - log 


^ 1^11 — ^10 — ^01 + ^00 1 ^ 


(S2) 


where Eab is the energy eigenvalue for \ah) and the distance between two excitations d is computed in the same 
way as in the previous section. We repeat this process for different logical ground states and disorder realizations; 
various values of Ug and A/t are shown as black squares (a-f) in Fig. S2(left). For each case, the histogram of inverse 
localization lengths is shown in the middle panel of Fig. S2. Note that we have chosen to histogram 1/^ rather than 
^ because the Jacobian induces long tails in 

The mean of the inverse localization length increases with both decreasing spin density and increasing disorder 
strength. It may seem that the variance of the distribution is also increasing with disorder strength in (a), (d), (e), 
and (f); however, this behavior does not continue at larger disorder strengths. Rather, as shown in the right panel of 
Fig. S2, the standard deviation saturates near A/t ^ 30 while the mean continues to increase. 


Initialization and DEER sequences — Here, we provide the detailed protocol for initializing the system into 
an effective many-body eigenstate using a combination of DEER (double electron electron resonance) spin echo and 
efficient quantum phase estimation [2, 3]. We assume that the quantum state of the system is initially in an unknown 
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superposition of many-body eigenstates. Our goal is to efficiently prepare a quantum state such that the coherence 
of qubit excitations and their interactions is maintained for sufficiently long times. In doing so, we only use local 
operators and projective measurements limited to the accessible region of size i. We assume that we can start by 
preparing the small local regions into their respective local eigenstates (e.g. eigenstates of spin 

operators within the region) [4]. 

The state of the full system can be written in the basis as 

IV^o) = |0 l) (8) |M) (g) |0i^) (S3) 

M 

where M enumerates all (exponentially many) possible configurations for the uninitialized regions. Then, we perform 
local spectroscopy in each region {L and R) in order to estimate the energy of a single excitation at the center of 
the region. Due to the diagonal interactions in the MBL phase, the spectral lines for these excitations are broadened 
by at most ^ texp(—//^ 2 )- With ^ ^ ^ 2 , one can resolve and manipulate individual localized excitations in each 
region within a timescale ^ We choose a single excitation in each region to use as our MBL qubits. These qubit 
excitations interact with spins in the rest of the chain and therefore dephase over a time scale T ~ t iexp(V6)- 
This dephasing can in principle be refocused by the spin echo sequence considered in Fig. S3a [2]. 


(a) Spin echo 



time 


(b) 


Left 


DEER 



TV 

Right " 

time 


FIG. S3. Spin echo pulse sequences, (a) The dephasing of the left MBL qubit can be refocused via a simple Hahn echo, (b) a 
modified spin echo sequence (DEER) isolates the dephasing of the left MBL qubit induced by interactions Sj = Slr + Slmr- 


For coherent interactions between MBL qubits, one needs to further initialize the system, since the interaction 
strength Sj depends on the configuration M of the uninitialized region. For a specific M, the interaction can be 
written diS dj = 6 lr ^ ^lmr [4]. Here, one can understand 6lr as the direct interaction between the two excitations 
while Slmr is the contribution provided by all of the multi-body interactions mediated by r-spins in the uninitialized 
region. It is this dependence of Sj on M that that prevents the coherent interaction between MBL qubits [2]. 
Interestingly, this problem can be solved by precisely measuring 6i (to within 6uj) since the measurement serves as a 
quantum non-demolition measurement for the configuration of the uninitialized region. Note that this simple solution 
is only possible due to the lack of thermalization in the MBL phase. 

For the measurement of (5/, we consider a modified spin echo sequence, known as DEER, as depicted in Fig. S3b 
[2]. This sequence refocuses the dephasing of the left MBL qubit induced by all but the MBL qubit on the right [2]. 
Each measurement of Sj projects the many-body state of the system into a superposition of a few eigenstates that 
share the values of Sj consistent with the measurement outcomes. Note that we do not require the preparation of 
a unique eigenstate as long as the prepared state enables coherent interactions between MBL qubits up to a desired 
time scale T ^ . Repeating the DEER sequence (followed by measurements of the left MBL qubit) with a free 

evolution time measures the phase shift modulo 27r. Therefore, by varying the free evolution time = 2~^ jduj 
with n = 0,1,2,..., one can efficiently measure dj up to the desired precision 5uj [3]. This particular sequence of 
Tn is chosen such that the protocol minimizes the total evolution time and is well described in Ref. [3]. The total 
measurement time ^Tn scales linearly with the inverse of the desired resolution Tinit = rjduj^ where a constant of 
order unity r] > 2 accounts for the repetition of the n = 0 step due to projection noise. 


Fidelity estimates — Here ,we analyze the fidelity of single qubit and gates under various conditions. The single 
qubit fidelity depends on both (i) whether the broadened spectral lines have power-law tails (e.g., a Lorentzian) or 
fall off more rapidly (e.g., a Gaussian) and (ii) how the spatial resolution I compares with ^i. When ^ < Ci, local 
level repulsion ensures that there are no levels within the local spectral gap £ of the targeted level that couple to the 
external driving field; when ^ > ^i, approximately 5uj/£ levels per localization length are within the spatial extent of 
the driving, where 5uj ^ I/T 2 is the spectral resolution. Moreover, when the linewidth is Lorentzian, the dominant 
imperfections are due to nearby transitions that use the tails ^ {5uj/£Y of the Lorentzian (spectral imperfection); 
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Gaussian 

Lorentzian 

^<6 

exp[—2£’/<5a;] 
(£/6) X {Slo/S) 

{Suj/Sf 

(£/6) X {5uj/8) 


TABLE L Fidelity of one-qubit gates for finite spectral and spatial resolution under various conditions. Two different types 
of spectral profiles arise from different origins: Gaussian from inhomogeneous broadening, including the effects of the bare 
un-enhanced interaction, and Lorentzian from homogeneous broadening, e.g., extrinsic decoherence processes. 


whereas, for rapidly decaying tails, the dominant imperfections are due to distant, exponentially weak lines that are 
within Scj of the targeted line (spatial imperfection). Fidelities in the four cases are summarized in Table I. 

We now focus on the fidelity of two qubit gates and derive Eq. (12) in the main text, which we present here again 
for convenience. 


where F is an effective decoherence rate and ^ 2 ^^ is an effective localization length for bus excitations. 
In general, the fidelity of a two qubit gate is given as 


1-^2 


max (r, ^/) 

^driven 


(S4) 


(S5) 


where the numerator is the rate of decoherence due to either external imperfections or bare dephasing interactions, 
while the denominator is the strength of the enhanced interaction. One can always optimize d such that ^/((iopt) = T. 
Then, using ^^nven ^ te ^opt/C 2 ^^ one can obtain Eq. (12). The dominant contributions to T and ^ 2 ^ vary under 
differing conditions. First, ^ 2 ^ is either ^2 or C ~ a^/tT 2 depending on the spectral resolution. When I/T 2 is 
larger than the typical energy spacing of bus-type excitations 82 , one necessarilly couples a single MBL qubit with 
multiple excitations. In this case, the spatial extent of the qubit is broadened in the same way as that of a bound 
state coupled to a multi-particle contiuum of effective mass m* ~ 1/ta^, leading to ~ ^T 2 /m* ~ a^JtT 2 [5-7]. 
The effective decoherence rate T is limited by the rate of population loss as well as the extrinsic decoherence scale, 
r = max (Floss, 1 /^ 2 ), where Tioss can be estimated using Fermi’s golden rule 


r 1/^2 

"" 1/T| + A2 



(S6) 


The first factor comes from the Lorentzian spectrum, the second from the effective coupling strength, and the last 
from the density of bus excitation states. In order to maximize the size of the dressed qubit, A needs to be as 
small as possible, limited only by I/T 2 . Therefore, we set F ^ I/(T 2 f 2 ) (V^ 2 )^, wherein the corresponding optimal 
spacing is given by dopt = In \tT 2 min (l, T 2 f 2 ^ 2 /^^)] • There is a cross over at P/T 2 = ^ 2^25 when /T 2 > £ 2^2 
loss error limits the fidelity, and dopt has to be relatively small in order to to avoid long gate operation times. For 
P/T 2 < £ 2^21 ^he spectral lines of the bus-excitations are not resolved, but the fidelity remains limited by the number 
of gate operations within a decoherence window. Thus, one finds 


1 -J-2 ^ (tT2)-l+«i/«'. 


(S7) 


We can estimate the fundamental limit of our protocol by considering a system optimized for our purposes viz. a system 
that supports different types of excitations with ^ 2/^1 ^ cx) with perfect spatial resolution and finite decoherence 
time T 2 < 00 . In this case the optimal fidelity is achieved from Eq. (S7). Consequently, in the limit ^2 ^ 00 and 

^ a, we obtain I — (tT 2 ) as presented in the main text. 
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